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We present improved measurements of the branching fractions of the color-suppressed decays 
B° _D'*'°/i'^, where h'^ represents a hght neutral meson tt", rj or lu. The measurements are 
based on a data sample of 140 fb~^ collected at the T(45') resonance with the Belle detector at 
the KEKB energy-asymmetric e"'"e~ collider, corresponding to seven times the luminosity of the 
previous Belle measurements. All the measured branching fractions fall in the range 1.4-2.4 xlO~'*, 
which is significantly higher than theoretical predictions based on naive factorization. 

PACS numbers: 13.25.Hw, 14.40.Nd 



I. INTRODUCTION 



The weak decays B° ^ D'-*^°h° where h° represents a light neutral meson, are expected to proceed predom- 
inantly through internal spectator diagrams, as illustrated in Fig. [T^. The color matching requirement between the 
quarks from the virtual W~ and the other quark pair results in these decays being "color-suppressed" relative to 
decays such as D^*^^h~ , which proceed through external spectator diagrams as shown in Fig.[T)3 . 



(a) (b) 




FIG. 1: Tree level internal (a) and external (b) spectator diagrams for B Dtt decays. 



Previous measurements of decays into Q and into D^*'>^Tr^, D^t] and D'^uj by the Belle collabora- 

tion, and of B° into D^*^°7r° by the CLEO collaboration [J, and into D^*^"tt^, D^*^"rj and D^*^"u; by the BaBar 
collaboration i aU indicate color-suppressed branchmg fractions in the approximate range (2-4) x 10 ^. Further 
color-suppressed branching fraction measurements of B^ decays into D^*^'^ri' by the Belle collaboration Q yield results 
of approximately (1.1-1.2) x 10~^. Most of these measurements are substantially in excess of theoretical expectations 
from "naive" factorization models [7-13] , which fall in the range (0.3-1.0) x lO^'*. 

Several approaches to achieving a better theoretical description [13, [l3, [IE] have been developed. They extend 
upon the factorization approach with consideration of final state interactions and consequent simultaneous treatment 
of isospin amplitudes of color-suppressed and color-allowed decays. The possibility that similar effects could have dra- 
matic implications on direct CP violation asymmetries in charmless decays, together with some degree of discrepancy 
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FIG. 2; Distributions of Mbc and AE for the decay modes B° D°tt° (a,b) and B° D*°tv° (c,d), with the three D° 
subdecay modes combined. The points represent the data, the solid lines show the result of the fit and the dash-dotted lines 
(peaking at AE « OGeV and Mbc ~ 5.28GeV/c^) represent the signal contributions. The vertical dotted lines represent the 
signal region. For the Mbc distributions the upper long dashed line shows the continuum-like background contribution, with 
peaking background contribution represented by the lower long dashed line and cross-feed contribution represented by the bold 
dotted line. For the AE distributions, (b,d) the two long dashed curves show the B background components (dominated by 
the color allowed B~ — > D'*^''p~ modes ) and the dotted line represents the continuum contribution and the bold dotted curve 
represents the cross-feed contribution. 



between the prior Belle Q and BaBar Q measurements provide strong motivation for more precise measurements of 
the color-suppressed decays. 

In this paper we report improved branching fraction measurements of decays into Z^^tt", D'^rj, D'^ui, £'*°7r°. D*^j] 
and D*^ijj. The measurements are based on a 140 fb~^ data sample, which contains 152 million BB pairs, collected 
with the Belle detector at the KEKB asymmetric-energy e+e" (3.5 on 8 GeV) collider [l^ operating at the T{AS) 
resonance. This corresponds to seven times the luminosity of the previous Belle measurements Q and almost twice 
that of the BaBar measurements Recent measurements of the same processes by the Belle Collaboration [l3| 
have been used to extract the angle 4>i of the CKM Unitarity Triangle using a time-dependent Dalitz analysis of 
D — > K'^TT^Tr~ , allowing resolution of the sign ambiguities inherent in other determinations. 

The Belle detector is a largc-solid-anglc magnetic spectrometer that consists of a three-layer silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold Cerenkov counters (ACC), a barrel- like 
arrangement of time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL) located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect mesons and to identify muons (KLM). The detector 
is described in detail elsewhere US . 
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FIG. 3: Distributions of Mbc and AE for the decay modes B° D°rj (a,b) and B° D'^rj (c,d), with the three D° 
subdecay modes combined. For the AE distributions, (b) the long dashed curve represents the sum of B background and 
continuum contributions; (d) the long dashed curve represents the B background contribution with the dotted line representing 
the continuum contribution; in both (b,d) the bold dotted curves represent the cross-feed contributions. In (b) two bold dotted 
curves are visible, corresponding to cross-feed arising from D*'^h° final states; the broader distribution is the component with 
final states where D*^^ and the other corresponds to final states where _D*° — > D'^tt''. Other conventions follow those 

described in the caption of Figure [2] 



II. EVENT SELECTION 



Color-suppressed i?" meson decays are reconstructed from candidate Z?" or mesons that are combined with 
light neutral meson candidates The 13" mesons arc reconstructed in the final states A'^7r+, A'^tt+tt". and 
A'~7r+7r"'"7r~ , while the light neutral mesons h'^ are reconstructed in the decay modes tt" — > 77, 77 —y 77. 77 tt^tt^tt^ 
and uj — > 7r+7r^7r°. The D*'^ mesons are reconstructed in the 13*° D°7r° decay mode. The decay mode D*^rj 
where ry — > Tr+Tr^Tr" is not reconstructed. 

Vertex and mass constrained fits arc performed for decays with charged products such as the three decay modes 
and 77 — *■ 7r"'"7r~7r'^; mass constrained fits are performed for the tt" 77 and rj candidates; and vertex constrained 

fits are performed for uj — *■ 7r"''7r~7r'' candidates. These kinematic fits result in reduced uncertainties on the energy 
and momenta of the candidate mesons. 

Charged tracks are required to have impact parameters within ±5 cm of the interaction point along the beam axis 
and within 1 cm in the transverse plane. Each track is identified as a kaon or pion according to a likelihood ratio 
derived from the responses of the TOF and ACC systems and energy loss measurements from the CDC. The likelihood 
ratio is required to exceed 0.6 for kaon candidates; within the momentum range of interest, this requirement is 88% 
efficient for kaons and has a misidentification rate for pious of 8.5%. 

The photon pairs that constitute 7r° candidates are required to have energies greater than 50 MeV and an invariant 
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FIG. 4: Distributions of Mbc and AE for the decay modes B° -> D°uj (a,b) and B" ^ D*°u (c,d) , with the three D° 
subdecay modes combined. The conventions follow those of Figure [S] 



mass within mass windows around the nominal ir^ mass ranging between ±2(t to ±3(t depending on the tt'^ momentum. 
The 7r° mass resolution (a) is momentum dependent, with values in the range a ~ 5.4—9.0 MeV/c^ over the momentum 
range of 7r° produced in the £)'^*^*'7r° decays considered. 

Candidate 77 mesons that decay to 77 are required to have photon energies i?^ greater than 100 MeV. In addition, 



the energy asymmetry 



1-^71 --^72 I 



is required to be less than 0.9. The 7] candidates are required to have invariant 

10.6MeV/c2 for the 77 ^ 77 mode and 3.4MeV/c2 
77 candidate are found to contribute 



-^71 +-^7 

masses within a 2.5(t mass window of the nominal mass, where a 
for the 77 — *■ 7r+7r^7r° mode. If either of the photons that constitute the 77 
to any 7r° 77 candidate, the 77 decay is excluded. The tt^ decay products of the 77 tt^tt'-k^ and uj —* tt~^tt~tt° 
candidates are required to have center of mass frame (CM) momentum greater than 200 and 500MeV/c, respectively. 
The Lu candidates are required to have invariant masses within zbST of the nominal mass value, where F is the natural 
width of the UJ meson (F = 8.49 ± O.OSMeV/c^ (H). 

Invariant masses of the candidates are required to be within ±2a of the nominal mass, where a is 8, 12 and 



and K 



5MeV/c2 for the K 
K-TT+ 

The D*^ meson candidates are obtained by combining candidate 



7r^7r° mode is required to be greater than 400MeV/c. 



modes, respectively. The CM momentum of the tt' in the 



and low momentum tt*^ mesons, where the 



soft 7r° momentum in the laboratory frame is required to be less than 0.6 GeV/c and the invariant mass difference 
\M{D*°) - M{D°)\ is required to be within 2MeV/c2 of the nominal value. 
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III. B RECONSTRUCTION 

The candidates are reconstructed from combinations of Z?" or D*'^ and /i" using the improved energy and 
momenta resulting from the vertex and mass constrained fits. 

Two kinematic variables are used to distinguish signal candidates from backgrounds: the beam-energy constrained 
mass Mbc = V(^boam)^ - IEp^P) and the energy difference AE = Y.E* - E*^^^, where E*^^^ is the CM energy, 
and E* , p* are the CM energy and momenta, respectively, which are summed over the or D*" and /i" meson 
decay candidates. 

The resolution of Aibc is approximately 3MeV/c^ for all modes, dominated by the beam energy spread, whereas 
the AE resolution varies substantially among modes, depending particularly on the number of tt'^'s in the final 
state. Candidates within the region \AE\ < 0.25 GeV and 5.2GeV/c^ < Mbc < 5.3GeV/c^ are selected for further 
consideration. Where more than one candidate or reconstruction hypothesis occurs in a single event, a consistency 
measure is used to pick the best candidate and reconstruction hypothesis. The consistency measure is constructed 
from the sum of the pcr degree of freedom for the relevant kinematic fits (to or and h^). For hypotheses 
including D*^, an additional term reflecting the deviation of the invariant mass difference \M{D*'^) — M{D'^)\ from 
the nominal value is included in the consistency measure. The fraction of events with multiple candidates with the 
same reconstruction hypothesis in the signal regions defined below is estimated from signal Monte Carlo samples to 
be less than three percent for all reconstructed modes. 

Signal region definitions in Mhc and AE are chosen based on their resolutions. A common A/bc signal region of 
5.27GeV/c^ < Afbc < 5.29GeV/c^ is used for all final states. The signal regions in AE' are mode dependent, with 
|A£:| < 0.08 GeV for D^t:^, D^tj{jj), D*°tt° and D**^!^ modes ; -0.08 GeV < AE < 0.05 GeV for D^rjin^mr) ; 
-0.12 GeV <AE< 0.08 GeV for D*°u; ; and \ AE\ < 0.05 GeV for D^w. 

The event yields and efficiencies presented in the following sections correspond to these signal regions. Figs. [21 [3] 
and m show the Mbc and AE distributions after application of all selection requirements and with the AE signal 
requirement applied for the Mbc distributions (a,c) and the signal requirement 5.27GeV/c^ < Mbc < 5.29GeV/c^ 
applied for the AE distributions (b,d). This selection includes continuum suppression requirements as described in 
the next section. The signal regions are indicated by vertical dashed lines on the figures. 

TABLE I: Measured signal region yields and MC estimates of signal region contributions for D'*'"/!" for the combined 

D" subdecay modes. The number of signal events (Neig) obtained from the AE fit are listed together with their statistical 
uncertainties. For the modes _D''*'°7r'' estimates of the background contributions from D^*^^ p (Nop) other B backgrounds 
(Nbbk), continuum (Nqq) and cross-feeds (A^^rs) are listed. Combined background estimates {N^kg) together with cross-feed 
contributions (Nxrs) are provided for the other D^/i" modes. For the other D*'^hP modes, background contributions are shown 
for B backgrounds (Nbtk), continuum (Nqq) and cross-feeds {N^crs)- 



Mode 


Ns^g 


Nbbk 


Nbkg 


^ xrs 




Ndp 




620.5 ± 39.1 


26.8 ±2.1 




25.9 ±6.5 


448.7 ± 28.0 


93.5 ±8.0 




160.7 ± 18.3 




131.8 ± 14.6 


6.9 ± 1.7 






^ 


64.7 ± 11.2 




51.2 ±8.2 


2.4 ±0.6 






225.6 ±21.5 
201.5 ± 20.1 




174.6 ± 16.0 
135.8 ± 12.9 


16.8 ±4.2 
12.8 ±3.2 






£,.0^0 


115.2 ± 14.9 


22.1 ±2.9 




2.0 ±0.3 


29.9 ± 3.9 


40.5 ±5.3 




49.8 ±10.0 


19.7 ±3.9 




2.2 ±0.4 


8.5 ± 1.7 






53.3 ±9.2 


26.4 ±4.5 




1.8 ±0.3 


19.5 ±3.4 





t The combined r] results are from a simultaneous fit to the individual r] samples, rather than summation of the individual yields. 



IV. CONTINUUM SUPPRESSION 

At energies close to the T(4S') resonance the production cross section of e'^e~ —^qq{q~ u, d, s, c) is approximately 
three times that of BB production, making continuum background suppression essential in all modes. The jet-like 
nature of the continuum events allows event shape variables to discriminate between them and the more-spherical 
BB events. 

Seven event-shape variables are combined into a single Fisher discriminant [l9j . These variables include the angle 
between the thrust axis of the B candidate and the thrust axis of the rest of the event (cos 9t), the sphericity variable, 
and five modified Fox- Wolfram moments. The technique and details of the variables used are provided in [19j. 
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Monte Carlo event samples of continuum qq events and signal events for each of the final states considered are 
used to construct probability density functions (PDFs) for the Fisher discriminant [l^ and cos^s, where is the 
angle between the B flight direction and the beam direction in the CM frame and additional angular variables in some 
modes, as indicated below. The products of the PDFs for these variables give signal and continuum likelihoods Cs 
and Cqq for each candidate, allowing a selection to be applied to the likelihood ratio C = Cg/ {Cs + Cqq). 

For the decays D*^ii° and D*'^rj, the vector-pseudoscalar nature of the decay products results in the longitu- 

dinal polarization of the I?*". The discrimination benefits from this polarization by incorporating PDFs for the 
helicity angle in the likelihood ratio. The D*'^ hclicity angle is defined as the angle between the direction of the 1?° 
and the opposite of the 5° direction in the D*'^ rest frame. The vector-vector nature of the decay products in the 
decay B^ D*^u! prevents the D*'^ helicity angle from being a useful discriminant in this mode, as the polarization 
of the decay products is not known. However, the uj tt^tt+tt" "splay" angle, defined as the angle between the 
directions of the 7r° and either the n'^ or n~ in the tt+tt" rest frame, is found to provide useful discrimination and 
is incorporated into the likelihood ratio. The method used to account for uncertainties arising from the unknown 
polarization is described in Section IVlIII 

Monte Carlo studies of the signal significance Ng/^/Ng + Nb, where Ng and Nb are Monte Carlo signal and back- 
ground yields (using signal branching fractions from previous measurements) , as a function of a cut on the likelihood 
ratio £ indicate a smooth behavior. Although the optimum significance is generally in the range 0.6-0.7, a looser cut 
of C > 0.5 is applied for all modes in order to reduce systematic uncertainties. This requirement removes (66-79)% 
of the continuum background samples while retaining (74-83)% of the signal samples. 

For the D^uj mode the polarized nature of the uj allows additional discrimination against backgrounds to 

be achieved with an additional requirement of |cos0;ie;| > 0.3, where the helicity angle 9hei is defined as the angle 
between the inflight direction in the w rest frame and the vector perpendicular to the to decay plane in the oj rest 
frame. 

TABLE II: Efficiency correction factors for the modes — > _D'*'"'/i°; the correction factors for combined modes are averaged 
using PDG subdecay fractions. 



Mode 






D^'iKwn'-') 


D"{Knnn) 




0.94 ±0.05 


0.98 ±0.05 


0.90 ±0.04 


0.98 ±0.06 




0.99 ±0.06 


1.03 ±0.06 


0.95 ± 0.05 


1.03 ±0.07 




0.94 ±0.06 


0.98 ± 0.05 


0.91 ± 0.05 


0.98 ±0.06 




0.97 ±0.06 








DOuj 


0.89 ±0.06 


0.92 ±0.06 


0.85 ± 0.06 


0.92 ±0.07 




0.86 ±0.09 


0.90 ±0.10 


0.83 ± 0.09 


0.90 ±0.10 




0.91 ±0.10 


0.95 ±0.10 


0.88 ± 0.10 


0.95 ±0.11 




0.83 ±0.10 


0.87 ±0.10 


0.80 ±0.09 


0.88 ±0.11 



TABLE III: Corrected efficiencies for the modes _D'*-"'/i'' excluding the subdecay branching fractions, for all modes 

combined and for the individual subdecay modes, as estimated for the AE fit samples. 



Mode 




D"{Kn) 




D"{Kmrn) 




0.075 ± 0.004 


0.173 ±0.008 


0.042 ± 0.002 


0.085 ±0.005 




0.061 ± 0.004 


0.140 ±0.008 


0.036 ± 0.002 


0.066 ± 0.004 




0.042 ± 0.003 


0.092 ± 0.005 


0.026 ± 0.002 


0.046 ± 0.003 




0.054 ± 0.003 


0.122 ±0.010 


0.032 ± 0.003 


0.058 ± 0.006 




0.025 ±0.002 


0.056 ± 0.004 


0.015 ±0.001 


0.030 ± 0.002 




0.036 ± 0.004 


0.088 ± 0.009 


0.019 ± 0.002 


0.042 ± 0.005 




0.039 ± 0.004 


0.093 ±0.010 


0.023 ± 0.003 


0.042 ± 0.005 




0.011 ±0.001 


0.026 ± 0.003 


0.007 ±0.001 


0.011 ±0.001 
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V. BACKGROUNDS FROM OTHER B DECAYS 

Significant background contributions arise both from color-favored decays and from other color-suppressed decays 
(cross- feed) D*^h^ . Some backgrounds have the same final state as the signal while others mimic signal due to 

missing or extra particles. 

Generic Monte Carlo [20| samples of BB and continuum qq are used to study the background contributions in 
the Mbc and AiJ distributions. The sample sizes correspond to approximately three times the expectations from the 
data sample analysed. The BB event sample excludes the color-suppressed modes under investigation. Signal mode 
samples for each of the decay chains considered are generated and reconstructed separately. They are used to estimate 
the cross-feed contributions between modes using the branching fractions measured here with an iterative procedure. 

The dominant cross-feed contributions to the D^h^ decays are found to arise from the corresponding D*^h^ decays. 
These contributions peak at the same A/bc as the signal but arc shifted to the lower side in /S.E. As can be seen 
from Figs, and [JJs, the cross-feed contribution is substantial in the region — 0.25GeV < AE < — O.lOGeV but 
quite small in the signal region. Cross-feed contributions to the D*^h^ decays are found to be small. In all cases, 
the fraction of cross-feed within the signal region is less than 10% of the observed yield. 

For the B^ Z3(*)°7r° modes, the dominant background is from the color-allowed B~ decay. In the 

B'^ D^TT^ mode non-reconstructed soft 7r° from Z?*° D°7r°, photons from D*'^ D'^j and tt^ from tt^tt*' 
produce the same final state as the signal. However, the missing particles cause a shift in AE with a broad peak 
centered at approximately AE = —0.2 GeV. In order to reduce contributions from this background, events that 
contain B candidates reconstructed as B^ —> D^*'>^p^ within the signal region 5.27 GeV < Afbc < 5.29 GeV and 
\AE\ < 0.1 GcV are rejected. This requirement reduces the color-allowed contribution in the region — 0.25GcV < 
AE < —0.10 GcV by about 60%; it docs little to reduce contributions in the signal region, but remains useful to 
facilitate background modelling. The A/bc distribution of these backgrounds is found to contribute at and slightly 
below the Afbc signal region; these are referred to as the "peaking background" . 

For the B^ decays to D^*^'^r] and D(*^°tt' modes there are potential backgrounds arising from non-resonant B^ 
TT+TT 7r° decays. Invariant mass, Af(7r+7r 7r°), distributions within the Afbc and AE signal regions indicate no 
significant contributions from these non-resonant decays. Ratios of data to Monte Carlo expectations in the invariant 
mass sidebands give values consistent with 1.0 with relative uncertainties in the range of 3-8 percent. As there are 
no indications of the presence of this background no systematic uncertainties from this source are assigned. 

TABLE IV: Measured branching fractions (xlO"") for the modes B° D'^'''>"h°, using separate subdecay mode samples, 
as obtained from the AE fit. 



Mode 










2.13 ±0.19 ±0.31 


2.02 ±0.24 ±0.33 


2.43 ±0.27 ±0.38 




1.91 ±0.29 ±0.24 


1.64 ±0.32 ±0.24 


1.66 ±0.40 ±0.23 




1.52 ±0.44 ±0.26 


1.51 ±0.46 ±0.28 


2.18 ±0.60 ±0.39 




2.57 ±0.37 ±0.27 


1.77 ±0.36 ±0.22 


2.60 ±0.42 ±0.31 




1.21 ±0.23 ±0.17 


1.45 ±0.36 ±0.23 


1.54 ±0.38 ±0.23 




1.43 ±0.42 ±0.25 


0.79 ±0.44 ±0.15 


2.09 ±0.57 ±0.38 




2.66 ±0.64 ±0.38 


1.90 ±0.62 ±0.31 


2.12 ±0.76 ±0.38 



VI. DATA MODELLING AND SIGNAL EXTRACTION 

Independent unbinned extended maximum likelihood fits to the AE and A/bc distributions arc performed to obtain 
the signal yields. The yields from the AE fits are used to extract the branching fractions, while the yields from 
the Afbc fits are used to cross-check the results. The principal parameters of the fits arc the normalization factors 
of the components used to model the observed distributions. The AE fit is performed in the range —0.25 GeV < 
AE < 0.25 GeV using the Afbc signal sample while the Afbc fit is performed in the range Afbc > 5.2 GeV/c2 using the 
mode-dependent AE signal samples. In most cases the shapes of the signal and background component distributions 
in Afbc and AE are obtained from fits to MC samples. 

The signal models used are the same for all modes, with the Afbc signals modeled with a Gaussian function and 
the AE signals modeled with an empirical formula known as the Crystal Ball (CB) line shape [2l|, that accounts for 
the asymmetric calorimeter energy response. The CB function is added to a Gaussian function in AE with the same 
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mean. Fits of the /S.E distributions to signal Monte Carlo for cacti final state are used to obtain the signal shape 
parameters. 

The cross-feed contributions in il/bc and A£' are studied using a combination of signal Monte Carlo samples from 
all other color suppressed modes, weighted according to the branching fractions obtained here. Smoothed histograms 
obtained from this combined sample are used as estimates of the cross-feed contributions. In the Aibc case the /S.E 
signal region requirement results in very small cross-feed contributions, which are fixed at the Monte Carlo expectation. 
For Ai? there are considerable contributions in the region —0.25 GeV < /S.E < —0.10 GeV. The normalization of this 
component is allowed to float in the fit for the ZJ^/i" modes; for the D*^h^ modes the cross-feed expectations are small 
and are fixed at the Monte Carlo expectation. 

Continuum-like backgrounds in the Mbc fits are modeled by an empirical threshold function known as the ARGUS 
function [2^ . The small peaking background contributions are modeled by a Gaussian of mean and width and 
normalization obtained by a fit to the BB background Monte Carlo Mbc distribution, using an ARGUS function plus 
a Gaussian. A systematic uncertainty of 50% is assigned to the determination of this small background distribution. 
This treatment allows the vast majority of the background to be simply modeled with the ARGUS shape, leaving a 
small but less well-known peaking background component that represents the deviation from the ARGUS shape. 

In fits to data il/bc distributions, the ARGUS background function parameters are fixed to the values obtained 
from fits to combined Monte Carlo BB and continuum background samples. The signal parameters are free, as are 
the normalizations of signal and background. The small peaking background and cross-feed contributions are fixed 
at their expected values. 

The /S.E background distributions in the B^ D^rj and B^ — > D^ui are modeled using smoothed histograms 
obtained from a combined continuum and generic BB Monte Carlo sample. For the i3° Z)(*)°7r*' modes, the shapes 
of the IS.E distribution arising from BB and continuum background are very different (see Fig. [J) necessitating 
separate modelling. The continuum shape is modeled with a first-order polynomial with slope obtained from fits to 
the continuum Monte Carlo sample. The shape of the BB background is modeled with a Gaussian function plus 
a second-order polynomial, with parameters determined from a fit to the generic BB Monte Carlo sample. In fits 
to data, the large peak in the region —0.25 GeV < Ai? < —0.10 GeV that arises principally from the color-allowed 
B~ — > D'^*^^ p" decays is found to be broader than the Monte Carlo expectation; thus all parameters of this color- 
allowed Gaussian are allowed to fioat in the fit. The normalizations of the contributions from the remainder of the 
BB background, the continuum and the signal are also floated in the fit, with the small cross-feed contribution fixed 
as discussed above. 

The results of the Mbc and /S.E fits for each of the D^*^^h^ modes with the three subdecay modes combined are 
presented in Figs.[2l [3] and [4] The D'-*^'^h^ mode results are obtained by a simultaneous fit to the three submodes, 
where the signal yields in each submode are constrained by the ratios of the products of efficiency and secondary 
branching fraction. 

TABLE V: Systematic uncertainties of the measured branching fractions for _B° _D'*''^ft°, for the combined D'^ submode 
samples, as estimated for the AE fit results. 



Category 


















Tracking efficiency 


2.6 


2.6 


2.6 


2.6 


2.6 


2.6 


2.6 


2.6 


/i" efficiency 


2.7 


4.0 


4.0 


4.0 


5.4 


2.7 


4.0 


5.4 


Kaon efficiency 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


Extra efficiency 


0.8 


0.8 


0.8 


0.8 


0.8 


0.8 


0.8 


0.8 


Likelihood ratio efficiency 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


MC statistics 


2.1 


2.1 


2.6 


2.3 


3.0 


3.3 


2.6 


5.3 


Slow 7r° (from D**^), efficiency 


0.0 


0.0 


0.0 


0.0 


0.0 


9.5 


9.5 


9.5 


Crossfeed 


1.1 


1.1 


1.1 


1.1 


1.6 


0.4 


1.1 


0.9 


Modelling ±lcr variations 


13.6 


11.1 


16.0 


8.3 


7.4 


12.9 


12.3 


6.4 


Subdecay Branching Fractions 


5.2 


5.2 


5.5 


5.3 


5.2 


7.1 


7.1 


7.1 


Number of BB events 


0.7 


0.7 


0.7 


0.7 


0.7 


0.7 


0.7 


0.7 


Longitudinal polarization fraction 
















7.1 


Total (%) 


15.6 


13.8 


18.1 


11.7 


11.9 


18.5 


18.3 


17.6 
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VII. BRANCHING FRACTION RESULTS 



The yields obtained from the A/bc and AE fits are consistent; the difference is typically within 50% of the statistical 
uncertainty. The results from the AE fits are found to have a slightly smaller total uncertainty in most cases and 
are used for the final result. The yields for the D^*^^hP modes (with the three D° subdecay samples combined) 
obtained from the one dimensional AE fits are presented in Table [D The statistical significances of the signals for 
each of the combined subdecay samples are greater than Qa. The combined t] ^ and rj — > 7r"'"7r~7r° yields are 
obtained from a combined fit to the individual 77 samples, rather than from a summation of the individual sample 
yields. The agreement between these approaches is apparent from the Table. For the D'^h'^ modes, backgrounds 
arising from cross-feed contributions in the signal region can be seen to contribute substantially less than the extent 
of the statistical uncertainty on the signal yield. 

The yields obtained are interpreted as branching fractions using the number of BB events, the product of subdecay 
fractions [2^ corresponding to the decay of D^h^ or D*^hP into the observed final states, and the efficiency. The 
efficiency for each mode is first obtained from signal Monte Carlo samples and then corrected to account for differences 
between data and MC expectations. 

The total efficiency corrections for each final state are obtained from the product of the relevant efficiency correction 
factors. The values for all final states of the D^*^'^hP modes are presented in Table HIl and the efficiencies are presented 
in Table Uni 

For the D*^h^ modes, the correction to the reconstruction efficiency of the soft pion produced in the process 
£)*o _^ IJOttO is estimated to be 0.920 ± 0.087; this dominates the correction value and uncertainty for these modes. 
This soft pion efficiency correction is estimated by comparing the yields of D*+ from the processes D*^ D^tt'^ and 
D*"*" D'^tt'^, using the subdecay modes 1?° K~7r~^, — > K~tt'^'k'^'k~ and D"*" K~'k'^tt^ . The final states 
differ by a single 7r° or 7r+, allowing the ratio of yields of D'^tt^ and D*+ D+tt" to be used to estimate 

the ratio of efficiencies between soft tt"*" and 7r°. Forming a double ratio of data over Monte Carlo and using separate 
estimates for the soft 7r+ efficiency correction provides the soft 7r° efficiency correction factor. The uncertainties on 
this factor are dominated by the uncertainties on the subdecay branching fractions. Reconstruction efficiencies for 
more energetic tt'^'s are obtained from comparisons of ry ^ tt^tt'^tt'^ to 77 — > 77 and to 77 ^ 7r"'"7r~7r°, for data and 
Monte Carlo. 



TABLE VL Measured branching fractions (xlO for the modes D^'^'^h'^, obtained from simultaneous fits to tlie three 

subdecay mode samples as obtained from the AE fit. 



Mode 


Branching fraction (xlO 




2.25 ±0.14 ±0.35 




1.77 ±0.20 ±0.24 




1.78 ±0.30 ±0.32 




1.77 ±0.16 ±0.21 




2.37 ±0.23 ±0.28 




1.39 ±0.18 ±0.26 




1.40 ±0.28 ±0.26 




2.29 ±0.39 ±0.40 



nO J-)Of,0 

TABLE VIL Ratios of branching fractions, £)*OfeO ' ^'^^^ statisticai and systematic uncertainties. 



Modes Ratio of branching fractions 

1.62 ±0.23 ±0.35 
^J-''^ 1.27 ±0.29 ±0.25 



"■^ 1.04 ±0.20 ±0.17 



The branching fraction results obtained from the AE fit yields for the D^*^'^h^ modes with the D° submodes 
combined are provided in Table IVTl Results from the individual submodes are shown in Table HVl These results are 
also graphically presented in Fig. O where the BaBar results from a one dimensional fit to Mbc distributions are 
also included for comparison. Table IVlTl presents the measured ratios of branching fractions, which benefit from a 
partial cancellation of systematic uncertainties. 
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TABLE VIII: Dependency coefficients expressing the relative clianges to tlie branching fractions provided in Table IVII due 
to relative changes to the branching fractions. The uncertainties presented are propagated from the uncertainties on the 
efficiencies only. 



Mode 


D'^{Kn) 








0.352 ± 0.034 


0.311 ±0.030 


0.337 ±0.032 


D°r)~,~, 


0.350 ±0.038 


0.329 ± 0.036 


0.321 ±0.035 




0.335 ± 0.037 


0.339 ± 0.038 


0.327 ±0.036 




0.344 ± 0.038 


0.333 ± 0.036 


0.323 ±0.035 




0.334 ± 0.044 


0.318 ±0.042 


0.347 ± 0.045 




0.368 ± 0.072 


0.285 ±0.056 


0.346 ±0.068 




0.362 ± 0.072 


0.319 ±0.064 


0.319 ±0.064 




0.356 ± 0.079 


0.341 ± 0.076 


0.303 ±0.067 



Table IVIIII provides the fractional dependence of the 5° branching fractions listed in Table IVII upon the 1)° 
branching fractions. These allow corrections to the measurements to be made to account for improved determinations 
of the branching fractions [H : B{D° K-tt+) = (3.89 ± 0.09)%, ^ A'-7r+7r+) = (13.9 ± 0.9)% and 

B{D^ /^"Tr+Tr+Tr") = (7.49 ± 0.31)%. The coefficients arc obtained from the subdecay branching fractions and the 
efficiencies for each subdecay. The fractional change to the branching fractions can be obtained from the product of 
the relative change of the subdecay branching fraction and the corresponding negated coefficient. 

VIII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties estimated for the and D*^h^ modes for the combined D'^ submode samples are 

provided in Table |Vl The systematic uncertainties are obtained by averaging those for the individual final states 
using the PDG subdecay fractions. The "Tracking efficiency" category accounts for charged products of 
only, estimated using an uncertainty of 1% per charged track. The efficiency" category corresponds to efficiency 
uncertainties for i:^ , -q and w. These uncertainties are obtained from the uncertainty on the efficiency corrections which 
is estimated from comparisons of ratios of various processes between data and Monte Carlo samples, as outlined in 
the previous section. The "Extra 7r°" category corresponds to the tt^ from the A'^tt+tt" process. 

The cross-feed uncertainty is estimated as 25% of the contribution from this source in the signal regions. Uncertain- 
ties arising from the background and signal modelling used are estimated from the changes in the yields as a result 
of ±1(7 variations on the model parameters. 

The "Longitudinal polarization fraction" uncertainty is relevant only for the vector-vector final state D*^u!. Al- 
though the analysis does not directly use angular information in these modes, the polarization can effect the orientation 
and hence momentum of the slow tt*^ from the D*^ resulting in changes to the total efficiency. Monte Carlo signal 
samples with various longitudinal polarization fractions (/l) are used to estimate the size of the effect. The results 
are quoted assuming = 0.5 ± 0.5; the variation in the efficiency over the full range of possible values is used to 
estimate the systematic uncertainty. 

The dominant contributions arise from the background and signal modelling, slow tt*' efficiency for D*^ modes, 
and the subdecay branching fractions. The total uncertainty is obtained by summing the individual uncertainties in 
quadrature. 

IX. CONCLUSION 

Improved measurements of the branching fractions of the color-suppressed decays D^tt", D'^r], D^u, D**^7r°, 

D*^r] and are presented. The results are consistent with the previous Belle measurements but with considerably 

improved precision due to the sevenfold increase in statistics. The individual results are consistent within two 
standard deviations with the BaBar measurements Q ; however it is notable that all the measurements are lower than 
those of BaBar. 

The measured values fall in the range (1.4-2.4) xlO"**, which is significantly higher than theoretical predictions 
based on naive factorization. This discrepancy indicates that either final state rescattering is significant, or else the 
assumption that the second Wilson coefficient 02 is real and process-independent is invalid 0, . 
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FIG. 5: Comparison of the measured branching fractions (xlO~*) of the modes D^*^'^h'^, for each of the three D'^ 

subdecay modes and for the combined mode samples, as obtained from the AE fit, with the Mbc signal region requirement 
applied. The shaded band indicates the combined submode result. The BaBar results from Afbc fits are also shown. 
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